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Abstract 

We have investigated the behaviour of 
the depairing current J e i p in ferromag- 
net /superconductor/ferromagnet (F/S/F) trilayers 
as function of the thickness d s of the supercon- 
ducting layers. Theoretically, Jd p depends on 
the superconducting order parameter or the pair 
density function, which is not homogeneous across 
the film due to the proximity effect. We use a 
proximity effect model with two parameters (prox- 
imity strength and interface transparency), which 
can also describe the dependence of the super- 
conducting transition temperature T c on d s . We 
compare the computations with the experimentally 
determined zero-field critical current J c q of small 
strips (typically 5 /jm wide) of Fe/Nb/Fe trilayers 
with varying thickness d^b of the Nb layer. Near 
T c the temperature dependence J c q(T) is in good 
agreement with the expected behaviour, which 
allows extrapolation to T = 0. Both the absolute 
values of J c o(0) and the dependence on d^b agree 
with the expectations for the depairing current. 
We conclude that Jd p is correctly determined, 
notwithstanding the fact that the strip width is 
larger than both the superconducting penetration 
depth and the superconducting coherence length, 
and that Jd P {d s ) is correctly described by the 
model. 
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I Introduction 

A still relatively little explored area of research 
in non-equilibrium superconductivity concerns 
phenomena involving spin-polarized quasiparticles. 
Pioneering work on spin-polarized tunneling in 
conventional s-wave superconductors was per- 
formed by Meservey and Tedrow |5J], who studied 
different ferromagnets (F) in F/AI2O3AI tunnel 
junctions and found that the tunnel current 
can show varying degrees of spin polarization. 
More recently, experiments were performed by 
different groups in order to establish whether 
superconductivity can be suppressed by injecting 
spin-polarized quasiparticles In these cases 

the combinations existed of a d-wave high-T c 
superconductor (XBa 2 Cu 3 07, with X = Y, Dy) 
and a fully spin-polarized ferromagnetic manganite 
(A . 6 7B .3 3 MnO3 with A = La, Nd and B = Ca, 
Sr), either with or without a barrier of a different 
oxide; measured was the change in the zero-field 
critical current density J c q of the superconducting 
films upon applying a current bias through the 
ferromagnet. The results are not fully conclusive, 
and certainly not quantitative. Although gener- 
ally a suppression of J c q was observed, heating 
effects could not always be fully ruled out since 
the manganites are highly resistive metals (see 
the discussion in ref. 0]), and the geometry 
did not always allow to determine the area of 
the current injection, and therefore the injected 
current density. Moreover, since J c q in high T c 
superconductors generally is not the depairing 



current Jdp but involves flux motion, J c q is not 
a direct measure for the amount of depression 
of the superconducting order parameter. Similar 
experiments have to our knowledge not been 
performed with combinations of conventional 
metals, although that would have some clear 
advantages. The interpretation of results would 
not be complicated by e.g. inhomogeneous cur- 
rents in the ferromagnet or anisotropic gaps in 
the superconductor; lithographic techniques could 
be brought to bear in order to have well defined 
superconducting bridges and injection contacts; 
and it should be possible to identify the effects of 
the spin-polarized quasiparticles on Jdp- 

Still, two points deserve special interest. The first 
is that, in planning such an F/I/S experiment, 
there is the potential problem of insufficient knowl- 
edge of the tunneling process. This was already 
apparent in the work of Meservey and Tedrow cited 
above Q since the experiments always showed 
a positive sign for the spin polarization, even in 
the cases of e.g. Co and Ni where a negative sign 
was expected. Recently, this was explained by 
demonstrating that the choice of barrier material 
can strongly influence and even reverse the spin 
polarization of the tunneling current with 
obvious consequences for the interpretation of the 
injection experiments. It may be advantageous 
to also contemplate an (F or N)/I/F/S config- 
uration; in this case the barrier is only used to 
increase the energy of the electrons coming from 
an N or F contact, while the polarization now 
takes place in a thin F layer between barrier and 
superconductor. The disadvantage here is that 
the F layer in connection with the superconductor 
will suppress the order parameter and therefore 
Jdp in the S layer. Still, since the proximity 
effect (which gives rise to the depression) for S/F 
systems is understood reasonably well 0], the 
effect on Jdp may also be quantifiable. The second 
point for consideration is that even in the case of 
conventional superconductors the determination 
of Jdp need not be straightforward. The difficulty 
lies in the fact that the superconducting bridge 
must have a width of no more than both the 
superconducting penetration depth A and the 
superconducting coherence length £. The first is 
needed to avoid current pile-up near the edges 
(as a consequence of screening of the self-field), 



the second is required in order to avoid vortex 
nucleation and flow, which gives rise to dissipation 
before Jd p is reached. These conditions can be met 
e.g. for Al, which has a Bardeen-Cooper-Schrieffcr 
(BCS) coherence length £ of about 1.5 fim, while 
A can also be made of the order of 1 /jm by making 
the film thin enough. For Al-bridges of less than 1 
fim wide it was shown by Romijn et al. that the 
measured Jdp agreed very well with the theoretical 
calculations by Kuprianov and Lukichev @] based 
on the Eilenberger equations and therefore valid 
in the whole temperature regime below T c . For 
a material such as Nb, with £o and A of the or- 
der of 50 nm, such agreement need not be expected. 

In this paper we show that, at least close to T c , the 
values of the zero-field critical current J c o measured 
on bridge-structured Nb samples are essentially the 
values expected for the depairing current. Further- 
more, we measure the depression of J c q in trilayers 
of Fe/Nb/Fe as function of the thickness dwb of the 
Nb layer. We compare the behaviour of J c o(dNb) 
with the behaviour of T c (<i/vfc), and also with cal- 
culations of the proximity effect and the pairbreak- 
ing velocity using a two-parameter formalism based 
on the Usadel equations. We find that J c o(dNb) is 
well described by the same two parameters which 
describe the behaviour of T c (dNb)- The conclusion 
is that the suppression of the depairing current as 
a consequence of the depression of the order pa- 
rameter in S/F structures can be well described by 
proximity effect theory, making (F,N)/I/F/S injec- 
tion experiments a distinct possibility. 

II Depairing current; theory 

Close to T c , the classical Ginzburg-Landau (GL) 
result for the temperature dependence of the de- 
pairing current of a thin film, under the assumption 
of a homogeneous superconducting order parame- 
ter over the film thickness, is given by : 

Jg L (t) = Jdf(0)(l-tf/ 2 , (1) 

with t — T/T c . The prefactor Jd p is of the order 
of H c /\, with H c the thermodynamic critical field, 
and will be given more precisely below. For arbi- 
trary temperatures, calculations were performed by 
Kupriyanov and Lukichev, who essentially solved 
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the Eilcnberger equations for a superconductor car- 
rying a current, with the velocity of the condensate 
leading to a phase gradient |8|. Their results re- 
cover the GL behaviour near T c : 
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Here, the constants have their usual meaning, N(Q) 
is the density of states at the Fermi level per 
spin direction, and x(Pd) is the G'orkov func- 
tion controlled by the 'dirt parameter' pd = 
(Hvp)/ (^TtkBTcle), with £ e the electronic mean free 
path. In the dirty limit, ( pd — > oo) x{pd) ~ * 
1.33^ e /£oj this becomes : 



J d G p L (0) = 1.26 [eN(0)v F A(0)} (3) 

which is equivalent to the expression given by 
Romijn et al. (TJ : 
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which can also be written in terms of experimental 
parameters as 



Jft (0) = 7.84 



{k B T c f 1 
? 2 Hvf(p£) p 



-,1/2 



(5) 



This way of writing also emphasizes the propor- 
tionality j9 L (0) oc \/l/p, since the product pi is a 



dp 

materials constant. At low temperatures the value 
of Jdp saturates, reaching a zero-temperature value 
of : 



J dp {0) = 1.491e7V(0)-^-A 3 / 2 (0) 

= 0.486 [eAT(0)u F A(0)] J , (6) 

V so 

with D = l/3vp£ e the diffusion constant. Com- 
parison with eq. ^ shows that the ratio between 
the saturation value and the GL-extrapolated value 
equals Jd P (0)/Jg L (0) = 0.385. 
In the case of F/S (or N/S) multilayers, the super- 
conducting order parameter is depressed near the 



interfaces, and this has to be taken into account 
in calculating Jdp- For this we use the proximity 
effect model, based on the Usadcl equations (dirty 
limit conditions), which was also used for calculat- 
ing the depression of T c with decreasing thickness 
of the superconductor ||. Details will be given 
in Appendix A but here we briefly introduce the 
main parameters of the theory. In principle, the 
shape of the order parameter on both sides of the 
interface depends on the bulk transition temper- 
ature T c o, on the coherence lengths £s,f, on the 
normal state resistivities ps,F, and on the trans- 
parency T of the interface. From the boundary 
conditions for the order parameter (see eqs. A4 ) 
it follows that, apart from T c0 , only two indepen- 
dent parameters are needed, the proximity strength 
parameter 7 and the transparency parameter 7;,. 
The value of 7 = (S,sPs)/{£,fPf) can be fully de- 
termined from experiment; the only free parameter 
is 76 (0 < 7b < 00), which is approximately con- 
nected to the transparency T (with < T < 1) 
by 
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As was shown in ref ||, in F/S systems, T can 
be quite low for a high magnetic moment in the 
F-layer, which is presumably due to the suppres- 
sion of Andreev reflections by the exchange split- 
ting of the spin subbands. Fig. |l| gives the results 
of some typical calculations, performed for the sys- 
tem Fe/Nb/Fe with the appropriate proximity ef- 
fect parameters 7 — 34.6 and 76 = 42 (see sec- 
tion IV). Shown is Jdp(t) for two different thick- 
nesses (ds/£s = 20,7.5), normalised on the bulk 
value Jdp lk {Q) as given by eq. |[ Note that this in- 
volves a factor {T c /T^ ulk f/ 2 . The thickness depen- 
dence of T c and the normalised depairing current at 
T = (see the insert of Fig. [I]) are quite different, 
with a much stronger depression of the depairing 
current at relatively high thickness of the super- 
conductor. This can be qualitatively understood 
by noting that T c is a measure for the maximum 
value of the superconducting order parameter in 
the layer, while the depairing current comes from 
an average over the layer thickness, which also in- 
volves lower values of the order parameter. 
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Figure 1: The temperature dependence of the 
normalised depairing current Jd P (t)/ J^p lk (0) of an 
F/S/F trilayer for S-laycr thicknesses ds/S,s — 
20 (upper), 7.5 (lower). Parameters typical for 
Fe/Nb/Fe were used, namely 7 = 34.6 and -ft = 42. 
Insert : thickness dependences of the normalised 
depairing current at T = 0, and of T c of the same 
trilayer. The black dots are at ds/£,s = 20, 7.5. 

Ill Experimental 

Samples were grown on Si(100) substrates, by DC 
sputtering in a system with a base pressure of 
10~ 9 mbar in an Ar-pressure of 6 x 10~ 3 mbar. 
Sputtering rates were of the order 0.1 nm/s for 
Nb and 0.03 nm/s for Fe. One series of samples 
consisted of trilayers Nb/Fe/Nb with Nb thickness 
rfjvb = 5 nm and the Fe thickness c?i? e varying be- 
tween 2 nm and 25 nm. These were used to de- 
termine the magnetization Mp e of the Fe layers 
in the presence of Fe/Nb interfaces with a com- 
mercial SQUID-based magnetometer. The behav- 
ior of Mp e versus d,Fe could be well described with 
a straight line, yielding a magnetic moment per Fe 
atom of 2.36 ps (mb being the Bohr magneton), 
slightly above the bulk value of 2.2 ps and a mag- 
netically dead layer per interface of 0.1 nm. This 
value is somewhat lower than reported for MBE- 
grown samples ||[l(|. Two other series of samples 
consisted of trilayers of Fe/Nb/Fe with dp e = 5 nm 
and varying d^b- One set was structured by Ar-ion 
etching into strips with a width w — 100 pm, the 
other into strips with a width w = 6 pm, or some- 
times 10 pm or 20 pm. In both cases the length 



between the voltage contacts was 1 mm. The first 
set (deposited in two different runs) was used for 
measuring T c {dNb), the second set for both T c {dNb) 
and J c (djsfb)- In all cases, the typical width of the 
resistive transitions to the superconducting state 
was 50 mK. 

Also measured were single films of Fe and Nb 
with different strip widths in order to establish val- 
ues for the specific resistivity pFe,Nb (at 10 K), for 
T c and for the upper critical field B C 2(T). On av- 
erage, we find pp e m 7.5 pftcm, pNb ~ 3.7 pflcm, 
T c = 9 K and S = -dB c2 / d T = 0.24 T/K, yield- 
ing Cgl(0) = y/<S> /{2i:ST c ) = 12.2 nm. This cor- 
responds to £s = 7.8 nm. No special precautions 
were taken to shield residual magnetic fields. The 
zero-field critical current I c was determined at dif- 
ferent temperatures T by measuring current (/)- 
voltage (V) characteristics. For this, a DC current 
was switched on for the time of the order of 1 s 
and the voltage recorded, to prevent heating via 
the contacts. All samples showed a clear transition 
from the superconducting to the normal state, with 
a large and almost instantaneous increase in voltage 
at I c . Upon detecting this rise, the current was also 
turned off since the sample then started to heat im- 
mediately. Most samples also showed a small rise in 
voltage prior to the major transition, probably due 
to vortex motion. We shall come back to this point 
in the discussion. Important for the theoretically 
expected value of Jd p (0) is the value of the resistiv- 
ity of the superconducting layer (see Eq. ||) . This 
value, pNb, was extracted from the normal state 
resistance R n at 10 K of the patterned samples by 
assuming that the Nb layer and the 10 nm thick 
Fe layer (pFe = 7.5 pftcm) contribute as parallel 
resistors. The resulting values for p^b are given in 
Table [l], together with the strip width w and T c . 
The values for the thinner films (around 50 nm) 
are somewhat larger than what we usually find for 
single Nb films, and approach that value for the 
thick films. 

IV Results 

Fig. H shows the measured values for T c (dNb) for 
both sample sets, with the the two types of open 
symbols denoting the two deposition runs for that 
set, and the solid symbols denoting the samples 
used for measuring J c q. The overall data spread 
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Table 1: Parameters of the Fe/Nb/Fe samples and 
the single Nb film used for the critical current mea- 
surements. Given are the thickness of the Nb layer 
d,Nb 1 the strip width w, the critical temperature T c , 
the calculated specific resistance of the Nb layer 
PNb and the Ginzburg-Landau extrapolated criti- 
cal current at zero temperature J^ L (0). 




Figure 2: T c of the different sets of Fe/Nb/Fe trilay- 
ers. The solid symbols denote the samples used for 
the critical current measurements. The line shows 
the theoretical dependence T c (d s ) for the parame- 
ter values 7 = 34.6 and 7b = 42. 
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Figure 3: Experimentally determined critical cur- 
rent density J c q versus reduced temperature t = 
T/T c for the Fe/Nb/Fe trilayers with d Nb = 42 nm 
(triangles), 60 nm (solid diamonds) and 75 nm 
(open squares). The solid and open symbols for 
dNb = 42 nm corresponds to measurements with 
non-pumped and pumped He-bath respectively. 



is small, and the data can be well described by the 
proximity effect theory for S/F systems we used for 
analyzing the behavior of V/(Fe x Vi_ x ) in ref. 0, 
with the two parameters 7 and 7b defined above. 
We use the same value for £p as in the case of 
V/Fe, £i?e = 0.14 nm and values for £ s , pp and ps 
as given in section H, yielding 7 = 34.6. The best 
description for T c (dNb) then is for 7b = 42, as shown 
by the drawn line in Fig. ||. The critical thickness 
for the S-layer for onset of superconductivity d^ r 
can be taken either from the lowest measured value 
for T c or from the extrapolated value of the calcu- 
lated curve, df r — 29 nm, corresponding to a ratio 
dfr/^s — 3.7, which is somewhat higher than in the 
case of V/Fe where we found 3.2. Apparently, the 
effect of ferromagnet on superconductor is slightly 
stronger in the Nb/Fe case, but this is not the issue 
of the current paper. 

In Fig. U J c0 = Ico/(wd) is plotted versus reduced 
temperature t = T/T c for djvb = 42, 60 and 75 nm. 
All curves show a clear upturn with decreasing tem- 
perature in the region close to T c , above t w 0.9. 
Plotting J c 2 /3 (£) versus t results in a straight line 
in this temperature regime, which can be extrap- 
olated to t = 0. The ensuing values for J^ L (0) 
are given in Table [l] for all samples, and comprise 
some of the main experimental results. They can 
also be used to normalize the data. Fig. ^ shows 
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Figure 4: ( J c0 / J^ L (0)) 2 / 3 versus t = T/T c for 
Fe/Nb/Fe trilayers with different thickness d s of 
the Nb-layer, as indicated. The drawn line indi- 
cates the GL-behavior, the dotted line is the result 
of the full calculation. 
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Figure 5: J% L {Q)p 1/2 of the Fe/Nb/Fe trilayers 
scaled on the value of the single Nb layer versus 
superconducting layer thickness d,Nb- The result of 
the model calculations for 7 = 34.6, 7b = 42 is also 
plotted (solid line) as well as the dependence of the 
critical temperature T c /T c q on d^b (dashed line) 
for the same parameters. 



[Jco(t)/Jtf) L (0)] 2/3 versus t together with the line 
1 — t (the GL-behavior) and the result of the full 
theoretical calculation, which is now independent 
of the parameters. All data collapse on the univer- 
sal curve above t = 0.9. At lower temperatures, 
the thinnest films (d s = 36, 40, 42. 53 nm) fol- 
low the full calculation quite closely, even down to 
t sa 0.6. The difference between the data of 36 nm 
and 40 nm is mainly due to the choice of the nor- 
malization value, and reflects the accuracy of that 
determination. For thicker films the first deviation 
progressively shifts to higher t. 

V Discussion 

The first point to be discussed is whether the mea- 
sured values of J c q agree with the theoretical esti- 
mates for Jd p . The absolute value of Jjj, L (0) can 
be calculated with Eq. ||. The materials constants 
for Nb are well documented 11 1; we use the values 
v F = 5.6 x 10 5 m/s and pi = 3.75 x 10~ 16 ftm 2 . 
Eq. I then yields for the Nb film Jg L Nb (0) = 
1.70 x 10 12 A/m 2 , which is quite close to the 
experimentally determined value of J^ L Nb (0j) = 
1.5 x 10 12 A/m 2 (see Tabic [t]). It is also in good 
correspondence with the data presented by Ando et 
al. |l^] on films with a thickness of 100 nm and dif- 
ferent strip widths between 0.1 /im and 10 pm, who 



found a fitted value J§ p L Nb {Q) = 1.26 x 10 12 A/m 2 . 
It appears that the depairing current is directly 
probed by the measurement of J c q. 
Next we consider the dependence of J% L {Q) on the 
superconducting Nb layer thickness d^b- As Eq. || 
shows, Jdp(0) is proportional to yY/pm,- Since 
PNb of the samples differs, this leads to some varia- 
tion in the expected value for Jd P (0) which can be 
taken into account by multiplying J ( § i (0) by p^ b - 
Normalizing this value to the single Nb film yields 
the dependence on d^b as shown in Fig. |. J% L (0) 
in the trilayers is clearly reduced with respect to 
the bulk Nb value and increases with increasing 
<i/V£>, but much more slowly than T c does. The cor- 
respondence with theory is good at low d^b, with 
some deviations above d^b ~ 75 nm. This coincides 
with the findings on the temperature dependence 
of Jco(t), shown in Fig. [| : for small djvb there are 
only small deviations in the whole measured tem- 
perature regime, for large dm the deviations are 
large below t = 0.9. This suggests that at high d^b 
the extrapolation for J^ L (0) leads to somewhat un- 
derestimated values. In essence, we conclude that 
the model used to describe the depression of T c in 
F/S/S trilayers also adequately describes the be- 
haviour of Jdp- 

This leaves one final point to be addressed. In or- 
der to determine the depairing current it is usually 
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understood that two conditions have to be fulfilled 
[Q,[l3) : the current has to be distributed uniformly 
over the strip, and the width w should be small 
enough to preclude vortex formation and motion. 
In terms of penetration depth X(t), strip thickness 
d s and Ginzburg-Landau coherence length £ (t) this 
means : 

W < X eff{t) = X(t), d s < A e //(i) 

= X 2 {t)/d s , d s > X eff {t) 
w < 4.4f (t) = 4.4£(0) / (1 - t) 1 / 2 . (8) 



The total current / including the edges is given by 



Estimating A(0) from A(0) = 1.05 x \Q- 3 ^/p /T c 
we find it ranges between 67 and 113 nm. Both 
conditions mean for all samples 1 — t < 10~ 4 , much 
smaller than the region where J c o(i) oc (1 — t) 3 ^ 2 , 
(Fig. ||) , and the question is valid whether the cur- 
rent is uniform, as has implicitly been assumed in 
the analysis. 

Qualitatively, current is expected to pile up at the 
edges of the strip in order to minimize the self-field 
inside. The edge current will then sooner reach 
the value of Id p - By using = Id p /(wd), this 
would lead to underestimating the real value of Jd P - 
From the close agreement between the experimen- 
tal and theoretical values this does not appear to be 
the case. Quantitatively, the situation can be as- 
sessed that the self-field of the sample is completely 
screened (B z — in the sample). The current dis- 
tribution is then given by |l4j ] : 



J(x) 



h 



ndy/W 2 - x 2 



W <x<W, (9) 



where It is the transport current through the sam- 
ple, x is in the direction of the width w of the film, 
x = in the middle of the film and 2W = w. Ac- 
cording to this formula, the current diverges at the 
edges of the film. It can be assumed, however, that 
the field penetrates over a distance d/2 from the 
edges, but is kept out of the rest of the sample by 
the screening current. Then, the current within 
d/2 from the edges can be set equal to Jd p L (0) and 
beyond d/2 it decreases according to Eq. Q. The 
following calculation can be done for the Nb film. 
The transport current It in the screened part of 
the strip can be calculated from Eq. ^ by using 



I/d = 2d/2jg L (0) + 



W-d/2 



w+d/2 irdvW 1 



= djg L (0) x 




\W\ 



W-d/2 

) 

-W+d/2 



(11) 



J(x = W-d/2) = jg L (0) 



(10) 



The ratio I/(wd) / J$ p L (Q), which can be calcu- 
lated from Eq. [0], gives the fraction of J^, L (0) 
which would be actually measured as the critical 
current under the given current distribution, where 
the depairing current is reached at the edges. It 
can be easily seen that it equals 1 when the cur- 
rent is uniform. For the Nb film with w and d as 
given m Table [TJ, Eq. [D] yields a fraction of 0.11, 
an order of magnitude below what is actually mea- 
sured. The conclusion is that J(x) is much more 
uniformly distributed than might be expected. The 
reason is probably that a field and moving vortices 
exist in the film, indicated by a voltage onset be- 
low the jump to the normal state. This breaks up 
the Meissner state, causes a much more uniform 
current distribution and allows the correct determi- 
nation of the depairing current over a much larger 
region than expected on the basis of the condition 
w < A,£. Still, the deviations of J c o{t) compared 
to the theoretical behaviour at higher dpjj, in the 
Fe/Nb/Fe trilayers may be due to the low value of 
A and inhomogeneities in the current distribution 
at these high thicknesses. At low thicknesses, there 
are two effects which increase A above the bulk 
value. Firstly, for d s < X (around dj^b ~ 75 nm) 
the effective penetration depth increases according 
to X e ff — X(0) 2 /d, and can become significantly 
higher than A(0). Secondly, the suppression of the 
order parameter as measured by the decrease of 
T c /T c0 results in a higher value for A(0). From that 
point of view the full agreement between the mea- 
sured and calculated values of J c o(t) at the lowest 
thicknesses is not surprising. 

VI Summary 

In this paper we have addressed the question of 
the value of the superconducting depairing current 
in F/S/F trilayers with varying d s , where the su- 
perconducting order parameter is inhomogeneously 
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suppressed by the pair breaking in the F-layers. 
The same model which is adequately describes the 
suppression of T c with decreasing d s with two pa- 
rameters (proximity strength 7 and interface trans- 
parency 76 or T) can also be used to compute 
the suppression of the depairing current. Measure- 
ments of the zero-field critical current J c o (as de- 
fined by the current where the resistance jumps to 
the normal state value) in thin strips of Fe/Nb/Fe 
show that the temperature dependence near T c is 
as expected for the depairing current. Also the ab- 
solute value of J C Q of single Nb films is close to the 
theoretically expected value and the measured sup- 
pression of J C Q in the trilayers follows the calculated 
behaviour. We conclude that the current distribu- 
tion is homogeneous and that the depairing cur- 
rent is measured, even though the strip widths are 
larger that the superconducting penetration depth 
and coherence length. Also, the proximity effect 
model correctly describes the shape of the order 
parameter, at least in the superconducting layer. 
These findings can be of use in experiments on the 
effect of injecting polarized quasiparticles. 
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Appendix A : calculation of J dp 

We assume that the dirty limit conditions are ful- 
filled in both S and F layers, so that the F/S bi- 
layer can be described by the Usadel equations. In 
the absence of a depairing current in the S layer, 
and in the regime of large exchange energy in the 
ferromagnet (E ex ^> &bT c ) these equations were 
discussed extensively by Buzdin et al. |l5|] (see also 
Demler et al. fl6|]). Here we rewrite these equation 
in fl-parametrization (F = sin 8, G = cos 9) and in- 
clude the pair-breaking effects by current along the 
S film : 



Zs-^Os(z) -co sin 6 s (z) + A(z) cos 6 s (z) = 0, 

(Al) 

d 2 

^ F —^e F {z)-ism6 F {z) = ^ (A2) 



A bx(T/T ) + irT r - sin S ) = 0, (A3) 

where u) n — 7r(2n + l)T/T c is the normalized Mat- 
subara frequency, uj — \u n \ + Q 2 cos 9(z), A is the 
pair potential in a superconductor normalized to 
ttT c , & = (HDs^ttTc) 1 / 2 , £ f = {hD F /2E ex y/ 2 
and D F ,s are the coherence lengths and the elec- 
tronic diffusion coefficients in F and S metals. 
Moreover, Q = £ s dx/dx is the normalized gradient- 
invariant superfluid velocity in the x-direction, with 
X the phase of the pair potential A. There are two 
sources of pair-breaking in the problem, the volume 
one by the current and the surface one by the fer- 
romagnet. The latter is described by the boundary 
conditions at the FS interface (z = 0) 

ts^-Os = iZf^-Of (A4) 
dz dz 

lbtF^e F = sm(9 s -e F ) (A5) 

dz 

where the parameter 7 = ps£,s / Pf£,f describes the 
strength of the suppression of superconductivity in 
S by the ferromagnet. 

The parameter "ft, describes the effect of bound- 
ary transparency (coupling strength) between the 
layers. In the NS case, when the decoupling is due 
to the presence of an additional potential barrier at 
the interface, jb = Rb/ PfCf, with Rb the normal 
state resistance of the N/S interface [Q. In the 
F/S bilayer there is no general microscopic deriva- 
tion for 76, combining the effect of exchange split- 
ting and an additional interface barrier. A sim- 
ple estimate is still possible, when the exchange 
splitting is the main cause for intransparency Q. 
Then 76 = (2/3)(Zf/frO < (1 - T A )/T A > where 
Ta is the transmission probability of scattering be- 
tween the majority and minority spin subbands, i.e. 
the probability of Andreev reflection. This process 
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is implicitly described by the boundary condition 
7h^F^^F = sin(#g — Op) since 8p is off-diagonal in 
spin indices. Here the brackets < ... > denote the 
Fermi surface averaging, which is generally propor- 
tional to the overlap area of the projections of dif- 
ferent spin subbands onto the contact plane (l^,|l9| . 
As a result, T4 drops roughly linearly (for spherical 
Fermi surfaces) as a function of E ex , both for bal- 
listic and diffusive interfaces po|| . The supercurrent 
density is given by 

J s (z,Q) = ^Q T y sin 2 S . (A6) 

Since the superconducting pair potential A and the 
Green's function S are suppressed by the super- 
flow Q, the dependence J S (Q) must be found self- 
consistently. In the well known spatially homoge- 
neous case [^lj the function J S {Q) behaves non- 
monotonously : the supercurrent J s increases with 
Q at small Q, then reaches a maximum and fi- 
nally drops to zero, when A is fully suppressed by 
current. The depairing current is defined as the 
maximum of J S {Q). A similar situation holds in 
the spatially inhomogeneous case considered here, 
with the difference that the solutions for 9(z) and 
A(z) of the proximity effect problem (eqs. A1-A3) 
should be calculated selfconsistently for a given Q 
using the boundary conditions at the FS interface 
(eqs. A4, A5). This problem is solved numerically 
by the method applied previously to NS bilayers 
and described in detail in p2fl . Then the local in- 
dependent supercurrent density J s {z,Q) is calcu- 
lated from eq.A6 by summing the solutions sin 2 9 S 
over u) n . Finally the density is averaged over film 
thickness J S (Q) = dj 1 J^ s J s (z,Q)dz and the de- 
pairing current is found from the maximum of the 
dependence of ( J s ) on Q. 
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